Relativistic spherical plasma waves 
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Tightly focused laser pulses as they diverge or converge in underdense plasma can generate wake 
waves, having local structures that are spherical waves. Here we report on theoretical study of 
relativistic spherical wake waves and their properties, including wave breaking. These waves may 
be suitable as particle injectors or as flying mirrors that both reflect and focus radiation, enabling 
unique X-ray sources and nonlinear QED phenomena. 
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The study of nonlinear oscillations in plasma 0, 3 is one of the great successes of plasma physics, both from 
understanding the underlying processes and from applications. Nonlinear waves, driven by intense lasers, offer a 
unique possibility to accelerate charged particles providing the longitudinal electric fields of the unprecedented 
strength (see [H and references cited therein). However the implementation of these schemes had to wait until the 
invention of Chirped Pulse Amplification which allowed for the generation of extremely short electromagnetic 
(EM) pulses, and until the understanding was gained of the mechanisms of injection of electron bunches into the 
plasma wave accelerating fields Q . Subsequently the process of electron acceleration by nonlinear plasma wake waves 
was demonstrated in a number of pioneering experiments Q- 

The nonlinear waves in plasmas can be used not only for the acceleration of charged particles, but they also can be 
used to generate ultra-bright X- and 7-ray radiation via the mechanism that is usually referred to as a Flying Mirror 
(FM) Q. This mechanism makes use of the fact that the nonlinear waves in the wavebreaking regime are characterized 
by the formation of "spikes" in the electron density, where the density tends to infinity. Such modulations of the 
electron density moving with the relativistic velocity and acting as a FM are able to reflect a counter-propagating EM 
wave in the form of high-frequency radiation Q . 

In this Letter we describe the basic properties of nonlinear oscillations in plasma, extending the analysis done by 
Dawson for the case of nonrelativistic spherical waves [H to the relativistic regime and discuss their applications. We 
notice that such waves demonstrate behavior typical for waves with continuos spectra, i.e., the phase mixing and 
consequent wave breaking, which can be used for electron acceleration and FM generation. As was shown in Ref. 
Q the breaking of a plasma wave provides means of injection of an electron bunch into the wakefield for subsequent 
acceleration. For spherical plasma waves, phase mixing is an intrinsic property and wavebreaking takes place at a 
specific distance from the center of the sphere, providing a way for controllable injection. The spherical wave traveling 
towards the center of the sphere in the wavebreaking regime is a good candidate for the FM. The spherical flying 
mirror (SFM) compresses and tightly focuses the incoming radiation, thus achieving the highest intensification factor 
among all the FM schemes. 

Assuming that ions are at rest, we utilize the model of a collisionless cold plasma, governed by Maxwell equations 
and by hydrodynamics equations of an electron fluid in the fixed ion background with homogeneous ion density no. 
We consider the case when all the variables that characterize the fields and the plasma depend only on the radius r 
and time t: 

d t n + r~ 2 d r (r 2 nv) = 0, (1) 
d t p + vd rP = -E, (2) 

v= P (i+p 2 r 1 ^, (3) 



-- 2 d r {r 2 E) = l-n, (4) 

where we normalized the time and space coordinate to and cuj~^, the electron density n to uq, the electron 
momentum p and velocity v to m e c and c, and the electric field E to m e u} pe c/ 'e, respectively, with uj pe = (47rnoe 2 /m e ) 1 / 2 
being the plasma frequency. 

We perform transformation from the Euler coordinates (r,i) to the Lagrange coordinates (ro,t) related to each 
other by r = rg + £(ro, i), where ro is the initial position of the electron fluid element and £(ro, t) is its displacement. 
Solving the continuity equation ((T|) we obtain for the electron density: 

n = n [J(ro,t)}~ , (5) 
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FIG. 1: (a) Phase plane [p, ro + f) for ro = 2 and h — 5, 10, 15, 20, 25. (b) Dependence of the oscillation period on the radius, 
Po = 5. 



where J(ro,t) = (1 + (?£/<9ro)(ro + £) 2 /r§ is the Jacobian of the transformation from the Lagrange to the Euler 
coordinates. The transformation has a singularity at the point where the Jacobian vanishes J(tq,£) — ► 0. This 
singularity occurs either when £ = — ro or for d£/ 'drQ = — 1, which is known as a condition of the wave breaking, when 
the electron density and the gradient of the electron velocity become infinite. 

The time dependence of the electron displacement, £ = £(ro,i), is determined from Eqs. @ - (21), which in the 
Lagrange coordinates take the form: 

d tP = -E, (6) 

d^=p(i+p 2 r 1/2 , (7) 

d t [(r +0 2 E} = (r +0 2 dtt (8) 
The solution of Eq. (J5J yields the expression for the electric field E 

E - 3(n, + 2 • ( ) 

Using this we find the integral of motion: 

6(r + 

where the constant h marks the trajectory in the phase plane (p, ro + £)■ Equation (fTU)) indicates that the oscillations 
can never reach the center of the sphere due to the singularity at £ = — ro. This also can be seen in Fig. la where we 
present the phase plane of electron motion, (p, ro + £)> which illustrates the nonlinear character of plasma oscillations. 
Using the integral of motion given by Eq. (| 10(1 we obtain the solution of Eq. ([7} in quadrature: 

u 

The period of oscillations can be determined from this equation by integrating over the entire oscillation cycle, i.e. 
by choosing the integration limits to be equal to £ m m and Cmaxi which are the minimum and maximum values of 
coordinate £ and determined from the p = condition in Eq. (QI 



£™„,™ax - 2 [2{h - 1) + r§/3] 1/2 cos [(6 T tt) /3] , (12) 

where = arctan |l — r^ 6 [2(/i — 1) + rp/3] 3 j and minus or plus corresponds to £ m j n or ^ max , respectively. For 
Po 2> 1 the period of oscillations T is 

T = 2(6 Po ) 1/2 , (13) 

where po is maximum amplitude of the momentum oscillation. In a weak nonlinearity limit we can solve the equations 
of motion. Assuming £ <C ro, and 9t£ < 1, we obtain from Eq. (|7J) 

5«£ + £ - f a /»-o - (4£ 3 /3^) + (3/2)£ (9 t £) 2 ■ (14) 
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FIG. 2: The dependences of density and longitudinal electric field on the radius in the case of the wake generation by the laser 
pulse. The field of the driving pulse is showed by a dashed line. 



It can be seen from the form of this equation that the nonlinear terms are of two types. The first type (the first two 
terms in the r.h.s) is due to the spherical geometry of the oscillations, while the second one (the third term) is due to 
the relativistic effects [Icj . 

This equation can be solved using perturbations theory by expanding in series the displacement £ and the oscillation 
frequency lu: £ = £o + £i + • ■ • an d ui = 1 + oj\ + .... The zeroth order solution to Eq. (fT4"|) is £o = « cos t with k being the 
oscillation amplitude, which in general depends on the coordinate tq. The frequency is determined from the condition 
that there are no secular resonant terms in the equations for £i, £2, We find that the oscillation frequency depends 
on the coordinate and amplitude as 



oj = 1 + K 2 /(l2rl) - 3k 2 / 16 - 



(15) 



The frequency grows while the coordinate ro decreases in accordance with the result of numerical calculation of 
oscillation period based on Eqs. (fTTj) and (fT2)) . which is plotted in Fig. lb. 

The frequency dependence on the coordinate means that the spherical plasma wave has a continuos spectrum. 
Considering the evolution of such waves we assume they have a form aexp(iip) with the eikonal ijj(ro,t). Then 
the wavenumber and frequency are expressed via the eikonal derivatives, k = d ra ip and to = —dtip- By virtue of 



crossdifferentiation we have dtk + d r 



dtd ro ip — d ro dtip = 0. This yields dtk — —d ro u. Since d ro uj does not depend 



on time we obtain k = fco — d TQ uj t. If the wave was initially converging, fco < and (duj/dro) < 0, or for initially 
diverging wave fco > and d ra ui > 0, then at the time, t = t turn , where 



tturn — kg/ d ro U>, 



(16) 



the wave will stop and then go in the opposite direction. The sign of the gradient of the nonlinear plasma wave 
oscillation frequency is determined by the 3D geometry factor or/and by the dependence of the wave amplitude k on 
the coordinate. 

The time of the wave breaking can be estimated following the results of Ref. [ll|. The displacement £ can be 
represented as a superposition of oscillations that are the harmonics of plasma frequency, which depends on the 
distance from the center of the sphere, rg: 



£i r o,t) = ^ £m(ro) exp [irm)(ro)t] 



(17) 



Differentiating £(ro,t) with respect to ro and substituting it into the wavebreaking condition 1 + d^/dr^ = 0, we 
obtain for the breaking time: 



t br = -lu/ (d ro Lu <9 t £) . 



(18) 
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FIG. 3: (Color on-line) Diverging and converging spherical wake wave: density profiles 

Assuming that for a converging plasma wave near the breaking point the frequency can be approximated as ~ 1 , 
then for dt^ — > 1 we have % r ~ To and for dt£, -C 1 the wavebreak time is tbr ~ J"o(c?f£) _1 . This means that in the 
relativistic case the first period of the wave breaks at tq ~ 1, i.e., at the distance of about a plasma wavelength from 
the center. This behavior of nonlinear plasma waves can provide a mechanism of controllable injection of electrons 
into a laser wakefield accelerator. 

Now we consider a spherical wake wave driven by an ultrashort EM pulse. We describe the pulse by introducing 
the ponderomotive force [J] in the r.h.s of equations of motion as f p = —d r ^f, where 7 = (1 + p 2 + a 2 ) 1 / 2 , a is the 
vector potential of the laser pulse and the pulse electric field has the form exp(— (7*0 — t) 2 )Erf(ro)/ro, and taking into 
account that d r = d ro /d rQ r. However if the plasma wave does not break inside the EM pulse and the oscillations are 
far from the center of the sphere we have d ro r « 1. The results of the numerical solution of the motion equations in 
the presence of the ponderomotive force of this form are presented in Fig. 2 for the converging pulse. Pronounced 
spikes in the electron density can be seen, meaning that the wake wave entered the wave breaking regime, which 
can be utilized in the relativistic SFM scheme. Since the 7-factor of a breaking plasma wave is determined by the a 
group velocity of the driving laser pulse (7^ rj w(a/2) 1 / 2 ) and the amplitude of the laser pulse increases as 1/r as it 
propagates towards the focus, it follows that 'Jsfm ~ 2~ 2 / 3 w 4 / 3 at ru ~ clq2~ 1 / 6 lj~ 2 / 3 . 

In what follows the results of 2D PIC simulations of intense tightly focused laser pulse interaction with underdense 
plasma are shown. The simulations were performed using the REMP (Relativistic ElectroMagnetic Particle) code 
[12j . The typical run utilized a simulation box with mesh grid size of A/ 20. 

In Fig. 3 the distributions of electron density illustrating the diverging (3a) and converging (3b) plasma wake waves 
are presented. In the diverging case a linearly polarized (out of the simulation plane) laser pulse with duration of 
lOA/c, introduced at the left boundary with a = 5 and the width of 20A was focused at x = 20A, resulting in a focal 
spot of about 3A. It interacted with a 150A x 160A hydrogen plasma with the density n/n cr = 1.56 x 10~ 2 . The laser 
pulse drives the diverging plasma wave. In the converging case a linearly polarized (in the simulation plane) laser 
pulse with a = 2, duration of 5A/c and width of 50A at the left boundary was focused at x = 150A. It interacted 
with a 85A x 110A hydrogen plasma with the density n/n cr = 5.6 x 10~ 3 . The converging plasma wave was produced, 
which breaks near the focus in accordance with the theoretical estimates. 

This breaking, converging wave acting as a SFM can not only up-shift the frequency but also can focus the reflected 
radiation tightly, further increasing its intensity. In order to illustrate this phenomenon computer simulations were 
carried out of a counterpropagating laser pulse interacting with such converging SFM. The linearly polarized (out of 
the simulation plane) counterpropagating laser pulse originated at the right boundary with a = 2, A/? = 4A, width of 
20 A and duration 5A/c. In Fig. 4 the z-component of the electric field is shown. The reflected pulse is focused into a 
tiny spot and its frequency is substantially increased. 

The typical intensity scaling is / ~ 7 3 /o when reflecting from a conventional (plane wave) FM Q, where / and Iq 
are the intensities of the reflected and incident pulses respectively. However in our case the situation with a SFM is 
different since it is able not only to reflect but also to tightly focus the incoming radiation into a focal volume that 
is 647 6 smaller than that of the incident wave. Assuming total reflection of the incident pulse the reflected pulse 
intensity scales as / = 2567 8 /o- This is the maximum possible intensification factor assuming a reflection coefficient 
R = 1. According to the results of Ref. the reflection coefficient R w (cj^/ujg) 2 /2j 3 at wave break, where uja and 
u) s are the frequencies of the driver and source laser pulses, respectively. If uj^ = u s then we have 

7= 1287 5 / . (19) 



This intensification factor exceeds all other obtained previously in different schemes using a FM (E, Il3|. 
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FIG. 4: (Color on-line) The reflection of a EM pulse from a converging wake wave. 

Highly efficient SFM can in principle enable the study of non-linear effects of Quantum Electrodynamics (QED) 
with present day laser systems, and SFM can possibly lead to ultra-bright sources of X- and 7-rays. For 7 = 15 
and 1 J energy of 1 /im wavelength laser pulse reflected by a SFM will give rise to 10 16 several keV photons. For 
higher intensities the energy of the emitted photons becomes of the order of the particle energy, so here the quantum 
recoil effects should be taken into account. Also at this point a process that has no classical analog appears: e + e~ 
pair production by a photon in a strong EM field. These non-linear QED processes paired together can lead to 
the electron-positron avalanche to occur [l4|. The parameters that characterize the probability of these processes 
Xe = \J {F^Pu) 2 /m e Es and x 7 = ' {F iiV hk v )' 1 / 'm e E s should be of the order of unity for the avalanche to start. Our 
estimates show that the pulse with intensity of Iq ~ 10 22 W/cm 2 should be intensified by a SFM with 7 = 10 for 
this effect to be seen. Here Eg = 1-32 x 10 16 V/m is the Schwinger field [HI, is the tensor of the EM field and 
p u , k v are the 4- momenta of the electron and photon respectively. The SFM scheme will also enable the study of 
the Schwinger process Jl5j. i.e., the production of e + e~ pairs from vacuum under the action of strong EM field. The 
estimates based on [ly] show that the intensity of 10 27 W/cm 2 of the reflected pulse will be enough to observe the 
pair production. Such intensity will require Iq ~ 10 20 W/cm 2 and SFM with 7 « 10. 

In conclusion, in this letter we studied the basic properties of the spherical oscillations in plasma. The produced 
spherical plasma waves demonstrate the behavior typical for waves with continuos spectra, i.e the phase mixing and 
consequent wave breaking. Since the breaking is due to the geometry of the oscillations, where and when it happens 
is determined by the properties of a plasma wave. Thus this wave can be used for controllable injection of electrons in 
the accelerating field of LWFA. In addition, the converging spherical wave may be suitable for a flying mirror, which 
both reflects and focuses incident radiation, enabling unique X-ray sources with applications to nonlinear QED. 
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